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This  report  summarizes  the  work  performed  at  the  Pennsylvania  State  University 
during  the  period  01/04/95  -  15/05/96.  The  Principal  Investigator  for  part  of  this  period, 
08/01/95  -  1 5/05/96  was  Dr.  Brian  L  Weiner,  while  the  Principal  Investigator  for  the 
other  part  of  this  period,  01/04/95  -  08/01/95  was  Dr.  Michael  Frenklach,  the  graduate 
student  Andrei  Kazakov  was  coinvestigator  during  all  of  this  period.  Two  studies  are 
included  in  this  report,  one  carried  out  by  Kazakov  and  Frenklach  and  the  other  by 
Weiner  and  Frenklach.  The  first  one  examined  the  relative  importance  of  acetylene  versus 
PAH's  as  major  species  influencing  soot  mass  growth  rates  according  to  the  kinetic  model 
developed  by  Frenklach  et  al.  The  other  is  a  preliminary  investigation  of  the  initial  steps 
of  a  reaction  pathway  by  which  C3H3  rsdicals  can  b©  conv6rt0d  by  acetylene  to 

PAH'S. 


SUMMARY  OF  ACCOMPLISHMENTS; 


Theses; 

A.  Kazakov  "A  computational  study  of  the  effect  of  pressure  on  soot  formation  in 
laminar  premixed  flames."  July  1996. 

Presentations; 

A.  Kazakov  and  M.  Frenklach,  'Modeling  of  Soot  Coagulation  and  Aggregation 
at  High-Pressure  Conditions',  Fall  Technical  Meeting  of  the  Eastern  States 
Section  of  the  Combustion  Institute,  Worcester,  MA,  1995,  pp.  415-418. 

A.  Kazakov  and  M.  Frenklach,  'Dynamic  Modeling  of  Soot  Particle  Coagulation 
and  Aggregation  in  High-Pressure  Laminar  Premixed  Flames',  Work-in-Progress 
Poster  Paper  140,  presented  at  26th  Symposium  (International)  on  Combustion, 
Naples,  Italy,  July  30,  1996. 
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SUMMARY  OF  RESEARCH; 


On  Roles  Of  PAH  And  Acetylene  As  Soot  Growth 

Species; 

A.  Kazakov  and  M.  Frenklach 


The  predictions  of  the  present  model  indicate  that,  for  both  1  and  10-bar  C2H4-air 
premixed  flames,  soot  mass  growth  rates  are  dominated  by  the  contribution  of  C2H2, 
while  the  PAH  condensation  on  soot  surface  is  negligible.  On  the  other  hand,  Howard 
and  co-workers  argued  in  their  recent  studies  (Marr  et  al.  1994,  Benish  et  al.  1995)  that 
soot  growth  proceeds  mainly  through  PAH  deposition  on  the  soot  surface.  Since  the  most 
recent  work  (Benish  et  al.  1995)  has  been  carried  out  for  atmospheric  ethylene-air  flames, 
it  is  of  interest  to  compare  the  present  model  predictions  with  the  analysis  of  the 

experimental  results  conducted  by  Benish  et  al. 

Benish  et  al.  studied  two  C2H4-air  laminar  premixed  flames  with  the  equivalence 

ratios  equal  2.1  and  2.4.  Tar  material  (PAH)  and  soot  were  collected  at  the  different 
locations  above  the  burner  and  weighted.  The  results  were  reported  in  mass  units  (mg/L) 
as  a  functions  of  distance  from  the  burner.  The  analysis  of  the  relative  contributions  of 
PAH  and  C2H2  to  soot  mass  growth  was  conducted  using  the  following  equations: 


dM. 


-  YC2H2-soof^,H2-soot^C2H2  +Ytar-soof^tar-soot"^ar 


dt 


(1.1) 


dM, 


tar 


dt 


-  Y C2H2-tar'^C2H2-tar'”C2H2  Y tar-soof^ar-soot^tar 


(1.2) 


where  Msoot  and  Mtar  are  the  masses  of  soot  and  tar,  respectively,  ©  collision  efficiency, 
Z  collision  rate  calculated  from  kinetic  theory,  ffic2H2  and  mtar  molecular  weights  of 
acetylene  and  tar,  respectively. 

The  assumptions  used  to  analyze  the  data  are  summarized  below: 
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1)  the  average  PAH  collision  diameter  is  estimated  to  be  10  A,  and  the  average  PAH 
molecular  weight  is  500  g/mol; 

2)  soot  mass  density  is  1.8  g/cm^,  all  particles  are  of  the  same  size;  the  particle  number 

density  is  equal  lO^^  cm-^  and  independent  of  the  height  from  the  burner; 

3)  the  temperature  is  equal  to  1650  K  and  independent  of  the  height  from  the  burner;; 

4)  the  acetylene  mole  fraction  is  equal  to  0.02  and  independent  of  the  height  from  the 
burner; 

5)  collision  efficiencies  of  C2H2-soot  and  C2H2-tar  are  equal. 


Upon  application  of  these  assumptions,  Eqs.  (1.1)  and  (1.2)  can  be  rewntten  as  (see 
Appendix) 


=  1.69x10  A^ootY  C2  H2  -soot  +1-0x10  Mjar A/sootY  tar-soot 


tYA/cogt,  _  gq  jy2/3  ,  .  1  A  ..  1  a8  jy 

dt 
dM, 


tar 


8  2/3 

=  7.8  X  10  A/jaj.Y 021^2 -tar  ~  1-0  x  10  tar-soot 


dt 


(1.3) 

(1.4) 


In  Eqs.  (1.3)  and  (1.4),  A^soot  and  A/tar  are  given  in  units  g/cm^,  time  is  in  seconds.  The 
rates  of  soot  and  tar  mass  groAvth  were  evaluated  by  the  differentiation  of  the 
experimental  curves. 

The  solution  of  the  system  (1.3)  -  (1.4)  is  obvious.  Performing  the  summation  of 
equations  and  using  the  assumption  (5),  one  can  obtain 


<^Afsnot  /<^^  +  <^A/tar/ 


YC2H2-tar  YC2H2-soot  7. g  x  10^  AY^a,  + 1.69 X  lO^AY^J, 


(1.5) 


and 


1.0  xlO  *  ('i/AYjoor  ,  /TA..  ia3  ,.r2l3.. 

Y  tar -soot  -  ^  ^2/3  V  dt  ^  ^soaiJ  C2H2-soot 


soot 


(1.6) 
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Using  the  described  analysis,  Benish  et  al.  concluded  that  tar  contributes  from  80  to  94  % 
of  soot  mass  increment  for  the  flames  studied. 

While  the  simulations  of  the  flames  considered  by  Benish  et  al.  were  not  carried  out 
due  to  the  lack  of  the  input  data,  the  modeling  results  for  a  similar,  1-bar  ethylene-air 
flame  studied  experimentally  by  Jander,  Wagner  and  co-workers  (Jander  1992)  are 
available  (Kazakov  et  al.  1995).  This  flame  has  the  same  stoichiometry,  /  =  2.1  and, 
therefore,  can  be  used  for  qualitative  analysis  of  the  observed  trends. 

According  to  the  original  experimental  studies  of  Howard  and  co-workers  (Marr  et  al. 
1994,  Benish  et  al.  1995),  the  collected  tar  material  included  the  PAHs  starting  from 
about  two  aromatic  rings  and  higher.  Therefore,  for  the  consistency  with  the 
experimental  analysis,  the  mass  growth  rate  of  tar  was  evaluated  using  the  computed  rates 
of  change  for  the  PAHs  starting  from  naphthalene.  The  computed  masses  of  “tar” 
material  and  soot  are  presented  in  Figure  la.  It  should  be  emphasized,  however,  that  the 
present  model  of  soot  mass  growth  considers  the  contributions  of  the  PAH  compounds 
starting  from  pyrene,  hence,  the  computed  soot  growth  rate,  in  fact,  did  not  have  any 
contribution  from  2-3  ring  PAHs,  and,  equally,  the  calculated  2-3  ring  PAH 
concentrations  were  unaffected  by  soot  formation.  Following  the  same  approach  as  it 
was  done  in  the  reduction  of  the  experimental  data,  the  rates  of  “tar”  and  soot  mass 
growth  are  obtained  using  the  computed  linear  flame  velocity  (Figure  lb).  Having  all  this 
information,  it  is  possible  now  to  evaluate  the  collision  efficiencies  defined  by  Eqs.  (1.5) 
and  (1.6).  These  results  are  presented  in  Figure  2a.  As  can  be  seen,  the  predicted  values 
agree  almost  quantitatively  with  the  data  reported  by  Benish  et  al.  Moreover,  the 
contributions  of  tar  and  C2H2  calculated  from  Eq.  (1.3)  (Figure  2b)  indicate  the 
“dominant”  role  of  tar  material  as  soot  growth  species,  in  agreement  with  the  conclusions 
made  by  Benish  et  al. 
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The  presented  numerical  test  demonstrates  that  the  model  of  soot  formation  which 
suggests  the  leading  role  of  C2H2  as  soot  growth  species  do  not  contradict  the 
experimental  results  of  Benish  et  al.  It  also  illustrates  the  fact  that  the  analysis  conducted 
by  Benish  et  al.  cannot  be  served  as  a  proof  of  the  dominance  of  tar  contribution  in  soot 
mass  growth. 

The  mechanism  of  mass  transfer  from  acetylene  to  soot,  expressed  by  Eqs.  (1.1)  and 
(1.2)  is  shown  schematically  in  Figure  3.  It  has  two  pathways:  the  “direct  route  via 
reaction  of  C2H2  with  soot  (I),  and  “indirect”  route  which  goes  through  tar  growth  (II- 
III).  As  can  be  seen  in  Figure  2b,  the  rate  of  tar  growth  is  substantially  lower  than  the 
rate  of  soot  growth.  As  was  discussed  above,  both  acetylene  and  tar  soot  growth 
hypotheses  can  explain  the  observed  trend.  On  one  hand,  the  direct  route  (I)  can  be 
dominant,  while  both  (II)  and  (III)  would  be  slow  to  explain  the  observed  low  rate  of  tar 
production.  On  the  other  hand,  if  the  direct  route  (I)  is  slow,  the  rates  of  (II)  and  (III)  are 
fast,  i.e.,  tar  is  formed  and  consumed  rapidly  resulting  in  low  net  tar  production.  Thus, 
the  information  available  from  the  experiment  appears  to  be  insufficient  to  prove  or 
disprove  either  point  of  view.  The  answer  depends  on  the  ratio  of  the  “channels”  (I)  and 
(II)  which  has  to  be  guessed  or  estimated.  The  conclusions  of  Benish  et  al.  in  essence 
follow  from  the  assumption  (5)  which  suggests  similar  reaction  probabilities  for  C2H2- 
soot  and  C2H2-tar  collisions.  The  analysis  of  the  modeling  results,  on  the  other  hand, 
indicates  that  the  kinetics  of  tar  growth  is  more  complex  than  that  described  by  Eq.  (1.2) 
and  cannot  be  treated  as  a  simple  C2H2-tar  reaction.  According  to  the  model  predictions, 
the  ultimate  tar  mass  comes  from  the  reactions  of  smaller  aromatic  compounds  with 
acetylene.  These  processes  are  strongly  reversible  which  causes  low  net  rate  of  tar 
production. 
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APPENDIX 


THE  DERIVATION  OF  MASS  GROWTH  RATES 


1.  Input  Data. 


T=\ 650  K  -  the  temperature; 

[C2H2]  =  1.47  A  10-7  mol/cm3; 
mtar  =  500  g/mol; 

"JC2H2  =  26  g/mol; 

[tar]  =  Mar/wtar  =  2  A  10-3  mol/cm^; 

c/tar  —  1 0  A; 

^/C2H2  =  3.4  A; 


®soot  ~  1'8  g/cmj; 


/^soot=  10*^  cm-3. 


^C2H2  - 


%RT 


=  1.16x10" 


cm/s; 


Ctar~  2.64x10  jjjjj/g. 


^soot 


6M 


V3 


soot 


^^Psoot'^soot' 


^  -ttA/  =TrV'*^ 


\2/3 


\^Psoot^ 


A^,=1.52xlO’J4-, 


4  ,.f2/3 
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2.  Rates. 


'"Cjh/c.h^-soo.  =  JCjH^]  =1.69x  io'a/S, 

^  5 

'",arZmr.»o.  =-^^oo.J»rM=  l.OOx 

+ ‘it«r^‘^CjH2-l«[C2H2][tar]  =7.8  X  10  M,;,, 
All  rates  are  given  in  units  g  cm”^  s~^,  the  units  of  Mtar  3itd  -A^soot  ^re  g  cm”^. 
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Figure  lb 
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Figure  2b 
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A  Possible  Pathway  For  PAH  Production  In 
Sooting  Flames  Involving  Neutral  Species.; 


B.  Weiner  and  M.  Frenklach 


Over  the  last  few  decades  there  has  been  much  debate  and  analysis  of  the 
mechanism  of  soot  formation  in  hydrocarbon  fuel  combustion.^  It  is  agreed  that  there 
are  various  stages  to  this  process,  starting  with  gaseous  precursor  formation 
followed  by  particle  nucleation,  coagulation,  agglomeration,  and  surface  growth 
and  oxidation.  However,  the  details  of  the  mechanism  have  not  been  definitively 
established  due  to  the  complexity  of  the  process  and  the  difficulty  of  the  process 
and  the  difficulty  in  obtaining  comprehensive  and  unambiguous  experimental 
data.  Both  ionic  and  neutral  species  have  been  observed  in  sooting  flames, 
however  under  normal  flame  conditions,  the  concentration  of  neutrals  is  seen  to 
be  much  higher  and  would  thus  seem  to  favor  a  free  radical  mechanism. 

The  objective  of  the  present  study  is  to  investigate  theoretically  a  model 
chemical  reaction,  representative  of  reaction  steps  that  may  play  a  role  in 
molecular  growth  of  Poly  Aromatic  Hydrocarbons  (PAH's)  in  hydrocarbon  flames. 
The  neutral  CsHa*  radical  is  found  in  high  concentrations  in  sooting  flames, 
thus  a  very  plausible  chain  of  reactions  that  can  serve  as  a  template  for  PAH 
growth  in  such  an  acetylenic  rich  atmosphere  is 

CjHj*  +  CjHj^CjHs* 
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As  the  formation  of  multi  ringed  aromatics  are  anticipated  to  follow  the  same  or  a 
very  similar  mechanism  to  the  formation  of  the  first  ring,  an  analysis  of  such  an 
initial  reaction  chain  will  be  a  guide  for  PAH  growth. 

The  calculations  were  carried  out  using  the  Gaussian  92  Quantum 
chemistry  program.  Molecular  geometry's  were  determined  by  finding  local 
minima  of  the  energy  with  respect  to  variations  of  the  nuclear  framework  of  the 
molecules  and  molecular  fragments.  The  energy  was  calculated  within  the 
Unrestricted  Hartree  Fock  (UHF)  approximation  using  a  6-31  g**  basis  set. 
Transition  states  between  various  isomers  were  determined  by  finding  saddle 
points  of  the  UHF  energy  functional  corresponding  to  one  negative  root  of  its 
Hessian  matrix.  At  these  stationary  geometry's  (local  minima  and  saddle  points) 
energy  corrections  were  calculated  using  spin  projected  Moeller  Plesset  Second 
Order  (PMP2)  perturbation  theory.  All  the  calculations  refer  to  energies  at  OK. 
Various  isomers  of  CgHj*  and  C7H7  •  rearrangement  pathways 

between  the  radicals  and  C7H7  •  studied.  The  most  energetically 

favorable  pathway  was  selected  that  leads  from  •  to  an  aromatic 
hydrocarbon  in  an  acetylenic  atmosphere.  The  initial  stages  of  this  pathway  are 
displayed  in  figures  1  and  2.  The  first  figure  shows  the  energetics  as  calculated 
within  the  UHF  approximation,  while  the  second  figure  the  energetics  that 
include  corrections  calculated  within  spin  projected  perturbation  theory,  PMP2. 

T-T  0 

The  energies  shown  are  with  respect  to  the  lowest  energy  isomer  of  ’ 
found  from  the  Gaussian  92  calculations.  The  first  relative  energy  shown  on  the 
far  left  of  figures  1  and  2  is  for  well  separated  •  and  ^2^2 ,  to  the  right  of 
this  we  display  the  isomer  of  •  and  its  relative  energy  that  is  formed  from 
the  initial  addition  reaction.  This  reaction  has  no  electronic  barrier  nor  does  the 
next  one,  that  involves  a  bond  formation  and  transport  of  the  radical  site. 
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producing  the  species  shown  third  from  the  right.  This  species  then  undergoes  a 
hydrogen  rearrangement  with  a  barrier  of  approximately  29.95  kcals,  as 
calculated  with  both  the  UHF  and  PMP2  methods,  to  form  the  radical  shown  on 
the  far  right  hand  side  of  figures  1  and  2.  The  last  stage  in  this  reaction 
sequence  is  the  only  one  with  a  barrier,  as  the  forward  activation  energy 
required  for  this  last  step  is  much  less  than  that  released  by  the  initial  steps  of 
this  reaction  sequence  there  is  sufficient  energy  available  for  this  reaction  to 
occur.  The  backward  activation  energy  barrier  can  be  seen  to  be  larger  than  the 
forward  one,  thus  one  expects  these  sequence  of  reactions  to  lead  to  the 
production  of  the  isomer  of  ^5^5  *  on  the  far  right  of  this  figure. 

Starting  from  this  isomer,  again  in  an  acetylenic  atmosphere  we  display  in 
figures  3  and  4  a  sequence  of  reactions  leading  to  an  aromatic  radical.  This 
time  the  energies  calculated  within  the  PMP2  approximation  are  somewhat 
different  from  those  calculated  within  the  UHF  one.  The  barriers  calculated 
within  the  UHF  ansatz  are  20.2,  24.3,  20  and  23.6  kcal/mol  above  the  energies 
of  the  starting  reactants  *^5^5  •  and  •  jh©  energy  of  20  kcal/mol  shown  for 
the  barrier  between  species  3  and  4  is  estimated,  as  a  fully  converged  transition 
state  has  not  been  obtained  in  this  case,  but  a  partially  converged  energy  of 
~20  kcal/mol  has  been  obtained.  The  barriers  calculated  with  the  PMP2 
correction  are  7.9, 13.3,  15  and  7.9  kcal/mol  above  the  energies  of  the  starting 
reactants.  The  barrier  between  species  3  and  4  is  again  estimated  this  time  as 
15  kcal/mol  in  line  with  the  genfal  decrease  of  the  PMP2  barriers  compared  to 
the  UHF  ones.  Note  that  for  the  transition  states  2  -3,  3-4,  4-5  and  5-6  the  actual 
activation  energies  are  larger  as  the  energies  of  these  stable  intermediates  are 
below  that  of  the  starting  species  as  shown.  The  largest  barrier  in  the  PMP2 
case  (which  is  the  more  reliable  as  it  is  a  more  accurate  approximation)  is  ~38 
kcal/mol,  though  this  might  be  an  overestimate  as  it  depends  on  15  kcal/mol 
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guessed  as  an  upper  limit  for  the  transition  state  3-4.  From  an  energy  point  of 
view  the  reaction  sequence  should  proceed  from  the  left  hand  side  species 
shown  in  figure  2  to  the  aromatic  radical  shown  on  the  right  hand  side  of  figures 
4,  as  approximately  81  kcals/mol  of  energy  is  released  in  the  first  sequence 
shown  (figures  2)  and  the  backward  activation  energies  are  greater  than  the 
forward  ones  in  both  sequences  (figures  2  and  4).  A  more  detailed  kinetic  study 
based  on  these  results  needs,  however,  to  be  carried  out  in  order  to  see  if  such 
a  conclusion  is  indeed  supported.  The  reaction  sequence  of  figures  3  and  4 
display  a  rather  novel  rearrangement  from  a  five  membered  (structures  2  and  3) 
to  a  six  membered  ring  (structure  5)  via  a  tricyclic  intermediate  (structure  4).  The 
bond  that  breaks  during  the  transition  state  4-5  is  marked  on  the  figures  with  a 
wavy  strike-through  symbol. 
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calculated  using  the  spin  projected  MP2  approximation.  The  energies  are  in 
kcal/mol  relative  to  the  energy  of  the  final  product  at  the  right  hand  side. 
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Figure  3:  Reactions  starting  from  C5H5  radical  plus  C2H2  followed  by  rearrangements 
calculated  using  the  UHF  approximation.  The  energies  are  in  kcal/mol  relative 
to  the  energy  of  reactants  at  stage  1 . 
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